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Abstract 

We present a new procedure to count the number of real zeros of a class of univariate 
Pfaffian functions of order 1. The procedure is based on the construction of Sturm sequences 
for these functions and relies on an oracle for sign determination. In the particular case of E- 
polynomials, we design an oracle-free effective algorithm solving this task within exponential 
complexity. In addition, we give an explicit upper bound for the absolute value of the real 
zeros of an A-polynomial. 
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1 Introduction 

Pfaffian functions, introduced by Khovanskii in the late ’70 (see m), are analytic functions that 
satisfy hrst order partial differential equation systems with polynomial coefficients. A funda¬ 
mental result proved by Khovanskii ([7]) states that a system of n equations given by Pfaffian 
functions in n variables defined on a domain U has finitely many non-degenerate solutions in 0, 
and this number can be bounded in terms of syntactic parameters associated to the system. 

From the algorithmic viewpoint, [5] presents a summary of quantitative and complexity 
results for Pfaffian equation systems essentially based on Khovanskii’s bound. The known elim¬ 
ination procedures in the Pfaffian structure rely on the use of an oracle (namely, a blackbox 
subroutine which always gives the right answer) to determine consistency for systems of equa¬ 
tions and inequalities given by Pfaffian functions. However, for some classes of Pfaffian functions 
the consistency problem is algorithmically decidable: for instance, an algorithm for the consis¬ 
tency problem of systems of the type /i(x) > 0,..., /^(x) > 0, fk+i{x) > 0,..., fi{x) >0, where 
X = {xi,... ,Xn), fiix) = and F) (I < i < /) and h are polynomials with integer 

coefficients, is given in m- This result allows the design of algorithms to solve classical related 
geometric problems (see, for example, M)- More generally, the decidability of the theory of 
the real exponential held (i.e. the theory of the structure Mexp = (1^;+) q “j Oj 1; exp, <)) was 
proved in [8] provided Shanuel’s conjecture is true. 

‘Partially supported by the following grants: PIP 099/11 CONICET and UBACYT 20020120100133 
(2013/2016). 
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In this paper, we design a symbolic procedure to count the exact number of zeros in a real 
interval of a univariate Pfaffian function of the type /(x) = F[x, (p{x)), where -F is a polynomial 
in Z[X,Y] and is a univariate Pfaffian function of order 1 (see O Definition 2.1]). The 
procedure is based on the construction of a family of Sturm sequences associated to the given 
function /(x), which is done by means of polynomial subresultant techniques (see, for instance, 
m)- As it is usual in the literature on the subject, we assume the existence of an oracle to 
determine the sign a Pfaffian function takes at a real algebraic number. Sturm sequences in the 
context of transcendental functions were first used in [TB] to extend the cylindrical decomposition 
technique to non-algebraic situations. In m. this approach was followed to count the number of 
real roots of exponential terms of the form p{x) + q{x)e^^^\ where p, q and r are real polynomials. 
Later in [9], the same technique is applied to treat the case of functions of the type F(x,e*), 
where F is an integer polynomial. 

A function of the form 

/(x) = F(x,e'^(")), 

where F and h are polynomials with real coefficients, is called an Fi-polynomial f|16ji. For 
these particular functions, we give an effective symbolic algorithm solving the zero-counting 
problem with no calls to oracles. To this end, we construct a subroutine to determine the sign 
of univariate Fi-polynomials at real algebraic numbers. Our algorithms only perform arithmetic 
operations and comparisons between rational numbers. In order to deal with real algebraic 
numbers, we represent them by means of their Thom encodings (see Section 12.2p . The main 
result of the paper is the following: 

Theorem 1 Let f{x) = F’(x,e^^^)) he an E-polynomial defined by polynomials F € 7^[X,Y] 
and h € Z[X] with degrees bounded by d and eoeffieients of absolute value at most H, and let 
I = [a, b] be a closed interval or I = M. There is an algorithm that computes the number of zeros 
of f in I within complexity (2dF [)'^°^^^. 

Finally, we prove an explicit upper bound for the absolute value of the real zeros of an 
Fi-polynomial in terms of the degrees and absolute values of the coefficients of the polyno¬ 
mials involved. This bound could be used to separate and approximate the real zeros of an 
Fi-polynomial. It provides an answer to the ‘problem of the last root’ for this type of functions. 
Previously, in [18], the existence of such a bound was established for general exponential terms, 
but even though it is given by an inductive argument with a computable number of iterations, 
the bound is not explicit. Algorithms for the computation of upper bounds for the real roots of 
functions of the type P{x,e^) or, more generally, P(x, trans(x)), with P an integer polynomial 
and trans(x) = e^, ln(x) or arctan(x) are given in |9| and [lOj respectively. 

The paper is organized as follows: in Section [J] we fix the notation and recall some basic 
theoretical and algorithmic results on univariate polynomials. Section [3] is devoted to the con¬ 
struction of Sturm sequences for the Pfaffian functions we deal with. In Section U we present 
our general procedure for zero counting. Finally, in Section [5l we describe the algorithms and 
prove our main results on Fi-polynomials. 
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2 Preliminaries 


2.1 Basic notation and resnlts 

Throughout the paper, we will deal with univariate and bivariate polynomials. For a polynomial 
F G Z[X, y], we write degx(T) and degy(F) for the degrees of F in the variables X and Y 
respectively, H{F) for its height, that is, the maximum of the absolute values of its coefficients 
in Z, and cont(T) G '^[X] for the gcd of the coefficients of T as a polynomial in Z[X][y]. 

Note that, if pi,P 2 G are polynomials with degrees bounded by di and d 2 , and heights 
bounded by Fli and H 2 , then H{pip 2 ) < (min{(ii, ^ 2 } + ^)HiH 2 . 

If / is a real univariate analytic function, we denote its derivative by f' and, for k > 1, its 
A:th successive derivative by 

For 7 = ( 70 ,... , 7 Ar) £ with 7 * 7 ^ 0 for every 0 < i < N, the number of variations 

in sign of 7 is the cardinality of the set {1 < f : 7 i-i 7 i < 0}. For a tuple 7 of arbitrary 
real numbers, the number of variations in sign of 7 is defined as the number of variations in 
sign of the tuple which is obtained from 7 by removing its zero coordinates. Given x G M and 
a sequence of univariate real functions f = (/o,... , Jn) defined at x, we write x(f, x) for the 
number of variations in sign of the {N + 1)—tuple {fo{x ),..., /Ar(x)). 

We recall some well-known bounds on the size of roots of univariate polynomials (see m 
Proposition 2.5.9 and Theorem 2.5.11]). 

Lemma 2 Let p = G C[X], au 7 ^ 0. Let r{p) := max{| 2 | : 2 G C, p{z) = 0}. Then: 

i) r{p) < 1 -|- max 

a) r{p) < (1+ ^ ^ 

\ o<i<d-i / 

We will also use the following lower bound for the separation of the roots of a univariate 
polynomial with integer coefficients (see m Theorem 2.7.2]); 

Lemma 3 Let p G Z[X] be a polynomial of degree d > 2, and ai, ■ ■ ■ ,ad be all the roots of p. 
Then 

min{|ai — : Oi 7 ^ aj} > d ~ {d + 1 )~^ H{p)^~‘^. 

A basic tool for our results is the well-known theory of subresultants for univariate polyno¬ 
mials with coefficients in a ring and its relation with polynomial remainder sequences (see [U 
Chapter 8 ]). 

Let F(X,Y) and G{X,Y) be polynomials in Z[X, T] of degrees d and e in the variable Y 
respectively. Assume e < d. Following (TJ Notation 8.33], for every — 1 < j < d, let SResj be 
the jth signed subresultant of F and G considered as polynomials in Z[A][y]. By the structure 
theorem for subresultants (see (TJ Theorem 8.34 and Proposition 8.40]), we have that 

SReSe-i = -Remainder 

where lc(G) is the leading coefficient of G and, for an index i with 1 < i < d such that SResj-i 
is non-zero of degree j: 
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• If SResj_i = 0, then SResj-i = gcd{F,G) up to a factor in Z[X]. 

• If SResj_i / 0 has degree k, 

Sjti-iSResk-i = —Remainder(sfctj_iSResj_i, SResj_i) 

and the quotient lies in Z[X][y]. Here, si denotes the Ith subresultant coefficient of F and 
G as defined in [TJ Notation 4.22] and ti is the leading coefficient of SRes;. 

We define a sequence of integers as follows: 

• no = d + 1, ni = d. 

• For i > 1, if SReSn;_i ^ 0, then nj+i = deg(SReSni-i). 

The polynomials 

Fi •— SRes^i^—i 

are proportional to the polynomials in the Euclidean remainder sequence 
G. Moreover, the following relations hold: 

— —1-^2+2 for i ^ 

where Gi G Z[X][y] for every i. 

2.2 Algorithms and complexity 

The algorithms we consider in this paper are described by arithmetic networks over Q (see 
[2]). The notion of complexity of an algorithm we consider is the number of operations and 
comparisons in Q. The objects we deal with are polynomials with coefficients in Q, which are 
represented by the array of all their coefficients in a pre-fixed order of their monomials. 

To estimate complexities we will use the following results (see m)- The product of two poly¬ 
nomials in Q[X] of degrees bounded by d can be done within complexity 0{M(d)), where M(d) = 
dlog{d) log log(d). Interpolation of a degree d polynomial in Q[A] requires 0{M{d) log(d)) arith¬ 
metic operations. We will use the Extended Euclidean Algorithm to compute the gcd of two 
polynomials in Q[X] of degrees bounded by d within complexity 0{M{d) log(d)). We will com¬ 
pute subresultants by means of matrix determinants, which enables us to control both the 
complexity and output size (an alternative method for the computation of subresultants, based 
on the Euclidean algorithm, can be found in PQ Algorithm 8.21]). For a matrix in its 

determinant can be obtained within complexity 0{n‘^), where ui < 2.376 (see [3l Chapter 12]). 

For a polynomial in Z[A], we will need to approximate its real roots by rational numbers and 
to isolate them in disjoint intervals of pre-fixed length with rational endpoints. There are several 
known algorithms achieving these tasks (see, for instance, |15j and the references therein). Here 
we use a classical approach via Sturm sequences. The complexity of the algorithm based on 
this approach is suboptimal. However, the complexity order of the procedures in which we use 
it as a subroutine would not change even if we replaced it with the one with the best known 
complexity bound. 


associated to F and 

( 1 ) 

1 ( 2 ) 
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Lemma 4 Let p € '^[X] be a polynomial of degree bounded by d and e € Q, e > 0. There is an 
algorithm which computes finitely many pairwise disjoint intervals Ij = {aj,bj] with aj,bj G Q 
and bj — Oj < e such that each Ij contains at least one real root of p and every real root of p lies 
in some Ij. The complexity of the algorithm is of order 0{d? \og{H{p)/e)). 

Proof. The algorithm works recursively. Starting with the interval J = {—(1 + H{p)), 1 + H{p)], 
which contains all the real roots of p (see Lemma [2]), at each intermediate step, finitely many 
intervals are considered. Given an interval J = (a, 6] with {p = 0} fl J / 0 and |J| > e, the 
procedure runs as follows: 

• Let c = and Jr = (c, b]. 

• If p{c) 7 ^ 0, let Ji = (a, c]. 

• If p{c) = 0 and c — e > a, let I = (c — e, c] and J; = (a, c — e]. If p{c) = 0 and c — e < a, 
take I = (a, c]. (Note that, in any case, I contains a real root of p and has length at most 

e.) 

• Determine, for each of the intervals Jr and J/, whether p has a real root in that interval 
or not. Keep the intervals that contain real roots of p. 

The recursion finishes when the length of all the intervals is at most e. The output consists of 
all the intervals of length at most e containing roots of p, including the intervals I appearing at 
intermediate steps. 

In order to determine whether p has a real root in a given interval, we use the Sturm sequence 
of p and p' (see [U Theorem 2.50]), which is computed within complexity 0(M(d) log(d)) by 
means of the Euclidean Algorithm. 

At each step of the recursion, we keep at most d intervals together with the number of 
variations in sign of the Sturm sequence evaluated at each of their endpoints. For each of these 
intervals, the procedure above requires at most 2d + 1 additional evaluations of polynomials of 
degrees at most d. Then, the complexity of each recursive step is of order O(d^). 

Since the length of the intervals at the kth step is at most , the number of steps is at 

most 1 + |~log( )]. Therefore, the overall complexity is 0{df’\og{H{p)/e)). □ 

In order to deal with real algebraic numbers in a symbolic way, we will use Thom en¬ 
codings. We recall here their definition and main properties (see [H Chapter 2]). Given 
p G M[Ai] and a real root a of p, the Thom encoding of a as a root of p is the sequence 
(sign(p'(a)),..., sign(a))), where we represent the sign with an element of the set {0,1,-1} 
Two different real roots of p have different Thom encodings. In addition, given the Thom en¬ 
codings of two different real roots ai and 0:2 of p, it is possible to decide which is the smallest 
between ai and a 2 (see (TJ Proposition 2.28]). 

For a polynomial p G M[A], we will denote 

Der(p) := (p,p',... ,p*''^®^^^) 


A useful tool to compute Thom encodings and manipulate real algebraic numbers is an effec¬ 
tive procedure for the determination of feasible sign conditions on real univariate polynomials. 
For pi,... ,ps G M[X], a feasible sign condition for pi, ... ,ps on a finite set Z C M is an s-tuple 
(ui, ..., as) G {=, >, <}^ such that {x € Z : pi{x)ai0,... ,Ps{x)as0} 7 ^ 0. 
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Lemma 5 (see /TH Corollary 2]) Given pQ,pi,... ,ps € M[X], pQ ^ 0, degp* < d for i = 
0 ,..., s, the feasible sign conditions for pi,... ,ps on {po = 0 } can be computed algorithmically 
within 0{sd^ log^{d)) operations. Moreover, if po has m roots in M, this can be done within 
0{smd\og{m) log^((i)) operations. The output of the algorithm is a list of s-tuples in {0,1, —1}®, 
where 0 stands for =, 1 for > and —1 for <. 

3 Sturm sequences and zero counting for PfafRan functions 

Following [3], we introduce the notion of a Sturm sequence for a continuous function in a real 
interval: 

Definition 6 Let fo : {a, b) ^ W be a continuous function of a single variable. A sequence 
of continuous functions f = (/o,..., /at) on (a, b) is said to be a Sturm sequence for /o in the 
interval (a, b) if the following conditions hold: 

If fo{y) = 0? there exists e > 0 such that fi{x) ^ 0 for every x € {y — e,y + e) C (a, 6 ), 
x^y, fo{x)fi{x) <0fory-e<x<yand fo{x)fi{x) >0ify<x<y + e. 

2. For every i = l,...,N-1, if fi{x) = 0 for x e {a,b), then /j_i(x)/i+i(x) < 0. 

3. fN{x) 7 ^ 0 for every x G (o, 6 ). 

Recalling that, for a given x G M, v{i,x) denotes the number of variations in sign of the 
{N + l)-tuple {fo{x ),..., /Ar(x)), we have the following analog of the classical Sturm theorem: 

Theorem 7 Theorem 2.1]) Let fo : {a,b) —>■ M 6 e a continuous function of a single variable. 
Let f = (/o,... j/at) be a Sturm sequence for fo in the interval {a,b) and let a < c < d < b. 
Then, the number of distinct real zeros of fo in the interval {c,d] is u(f, c) — u(f, d). 

The aim of this section is to build Sturm sequences for a particular class of PfafRan functions 
we introduce below. For the definition of PfafRan functions in full generality and the basic 
properties of these functions see, for instance, [3]. 

Given a polynomial G Z[X,Y] with degy(d>) > 0, let p he a, function satisfying the 
differential equation 

ip'{x) = ^{x,ip{x)). (3) 

Note that p is analytic on its domain, which may be a proper subset of M. 

We are going to work with PfafRan functions of the type 

f{x) = F{x,ip{x)), 


where F G Z[X,T], 

Taking into account that the first derivative of such a function is 

OF dF 

— {x,p{x)) + —{x,p{x)).^{x,p{x)), 


we define, for any F G Z[X, T], the polynomial F G Z[X,Y] (associated with $) as follows: 
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Thus, we have that 


fix) = F{x,ip{x)). 

Due to the following result, in order to count the number of real zeros of a function / as 
above, we will assume from now on, without loss of generality, that Resy(T, T) ^ 0. 

Lemma 8 Let be as in equation (Oj and let F € Z[X,Vj with degy(F) > 0. There exists 
a polynomial P € Z[X,Vj such that Resy(P, P) / 0 and P{x,ip{x)) has the same real zeros as 
F{x,ip{x)). Moreover, the polynomial P can be effectively computed from F and 

Proof. Without loss of generality, we may assume that F is square-free. Suppose that Resy (P, P) = 
0. Write P = cont(P)Po. Then, Resy(Po,Po) = 0 and so, the greatest common divisor of Pq 
and Pq is a polynomial S G Z[X,Y] of positive degree in T. If 

Fo = SU and Fq = SV 


for U,V £ Z[X, y], we have that 

fo{x) = Fo{x,ip{x)) = S{x,ip{x))U{x,ip{x)) and /o(x) = Fo{x,ip{x)) = S{x,ip{x))V{x,ip{x)), 

which implies that a zero ^ of /o which is not a zero of U{x,ip{x)) satisfies that mult(^,/o) = 
mult(^, 5(x, y?(x))) < mult(^,/Q), leading to a contradiction. Then, /o and U{x,ip{x)) have the 
same zero set in M. As 

Fo = (PP) = SU + SU, 

it follows that, if P G Z[X,Y] is a common factor of U and U with positive degree in Y, then 
T divides Fq = SV. Since U and V are relatively prime polynomials, then T divides S and, 
therefore divides Pq, contradicting the fact that Pq is square-free. 

The lemma follows considering the polynomial P = cont(P) U. □ 

We will apply the theory of subresultants introduced in Section [2] in order to get Sturm 
sequences for /. 

Let ^ 

Pi = Remainder(lc(P)^P,P) G ^[X][Y], 

where D is the smallest even integer greater than or equal to 1 -|- degy (P) — degy (P). 

Notation 9 Following Section \2.1[ for i = 0,..., N, let Ri := SRes^^-i G Z[A][y] be the 
{ui — l)th subresultant polynomial associated to F and Pi, Ti := Pi-i S the leading 

coefficient of Ri and, for i = 2,..., N+ 1, let pi := Sm £Z[X] be the Uith subresultant coefficient 
of F and Pi. 

Definition 10 For an interval I = (a, b) containing no root of the polynomials Ti for i = 
0,..., N or Pi for i = 2,..., N + 1, we define inductively a sequence {o'i^i)o<i<N £ {1, —1}"^'''^ 
as follows: 

• <7/,0 = <7/,l = 1 , 

• Cj /2 = ( —l)l(^®Sv(P-degy(Fi))(degy(F)-degy(Fi)+l)g^^|-|^|-^^^^degy(F)-degy(Fi)+l^ 
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• 0-/,i+2 = Sgi{pi+2Ti+lPi+l'Ti)(y I ,i, 

where, for a continuous function g of a single variable with no zeros in I, sgj{g) denotes the 
(constant) sign of g in I. For i = 0,..., N, we define 

Fi,i = o-piRi eZ[X,Y]. 

Finally, if I is contained in the domain of p, we introduce the sequence of Pfaffian functions 
f/ = {fi,i)o<i<N defined by 

fpi{x) = Fpi{x,p{x)). 

Proposition 11 Let F € ZlX,Y], degy(F) > 0, and let p be a Pfaffian function satisfying 
p'{x) = where € 'L\X,Y\ with degy(<i>) > 0. Consider the function f{x) = 

F{x,p{x)). Let F G Z[X,Y] be defined as in &■ Assume that the resultant Resy(F, F) G 
Z[X] is not zero. With the notation and assumptions of Definition[TU. the sequence of Pfaffian 
functions f/ = {fi,i)o<i<N is a Sturm sequence for f in L = (a, b). 

Proof. For simplicity, as the interval L is fixed, the subindex / will be omitted throughout the 
proof. 

First we prove that /o and fi do not have common zeros in L. Suppose a G / is a common zero 
of /o and fi. Then F{a,p{a)) = 0 and Fi{a,p{a)) = 0; therefore, pN+iic() = Resy (F, Fi)(a) = 
0, contradicting the assumptions on I. 

From this fact, taking into account that fo = f, and fi has the same sign as f' at any zero 
of / lying in J, condition 1 of Definition [U] follows. 

To prove that condition 2 holds, first note that if fj{a) = 0 and /j+i(a) = 0 for some a € L, 
since pi and Tj do not have zeros in /, by identities ([T]) and ([2]), a is a common zero of all fs, 
contradicting the fact that /o and fi do not have common zeros in L. Then, condition 2 in 
Definition [ 6 ] follows from the definition of the signs Uj and identities ([T|) and ([ 2 ]) . 

Condition 3 follows from the assumption that ttv, which equals /tv up to a sign, does not 
have zeros in /. □ 


In order to count the number of zeros of a Pfaffian function in an open interval, provided 
that the function is defined in its endpoints, we introduce the following: 


Notation 12 Let / : J ^ M be a non-zero analytic function defined in an open interval J C M 
and let c G J. We denote 


and 


sg(/, c+) 


sign(/(c)) if /(c) / 0 
sign(/(^)(c)) if mult(c,/) = r 


sg(/,c ) 


sign(/(c)) if /(c) / 0 

sign((—l)’'/(^)(c)) if mult(c,/) = r 


where mult(c,/) is the multiplicity of c as a zero of f. 

For a sequence of non-zero analytic functions f = (/o,..., /tv) defined in J , we write c(f, c^) 
for the number of variations in sign in (sg(/o, c''), ... , sg(/Tv, c'')) and c(f, c“) for the number 
of variations in sign in (sg(/o, c“), ..., sg(/Tv, c“)). 


Note that sg(/, c"*") is the sign that / takes in (c, c + e) and sg(/, c ) is the sign that / takes 
in (c — e, c) for a sufficiently small e > 0. Then, by Theorem [3 we have: 

Proposition 13 With the assumptions and notation of Proposition \lll if, in addition, the 
closed interval [a, 6 ] is contained in the domain of ip, the number of zeros of the function f in 
the open interval I = (o, 6 ) equals — n(f/, 6 “). 

As a consequence, we get a formula for the number of zeros of the Pfaffian function / in any 
bounded interval: 

Theorem 14 Let f{x) = F{x,ip{x)), where F G Z[X, T], degy(-F) > 0, and ip is a Pfaffian 
function satisfying ip'{x) = <f>(x, (/?(x)) for a polynomial € Z[X,Y] with degy(<f>) > 0. Assume 
Resy(F, F) 7 ^ 0. Consider a bounded open interval (q;,/3) C M such that {a,j3] is contained in 
the domain of (p. 

Let Pi and Ti he the polynomials in Z[A] introduced in Notationl^ If ai < 02 <■■■< ak 
are all the roots in (a, /3) of pi and Ti, the number of zeros of f in [a, /3] equals 

k 

#{0 <j<k + l: f{aj) = 0} + ^ ^(f/. ,aj) - u(f/., 

j=0 

where oq = ol, Ofc+i = /3 and, for every < j < k, Ij = {aj,aj+i) and ij. is the sequence of 
functions introduced in Definition{T^ 

4 Algorithmic approach 

Let y? be a Pfaffian function satisfying 

ip'(x) = <I>{x,p>{x)) 

for a polynomial 4) G Z[A, Y], Let 5y := degy(<l)) > 0 and 6x ■= degjf (<1>). 

In this section, we describe an algorithm for counting the number of zeros in a bounded 
interval contained in the domain of of a function of the type 

fix) = Fix,ip{x)), 


where F G Z[A, Y] with degy(F) > 0. 

To estimate the complexity of the algorithm, we need an upper bound for the multiplicity of 
a zero of a function of this type. Here, we present a bound in our particular setting which takes 
into account the degrees in each of the variables X and Y of the polynomials involved in the 
definition of the functions. A general upper bound on the multiplicity of Pfaffian intersections 
depending on the total degrees of the polynomials can be found in [5l Theorem 4.3]. Even 
though both bounds are of the same order, our bound may be smaller when the total degrees 
are greater than the degrees with respect to each variable. 

Lemma 15 With the previous notation, let gix) = G(x,(/?(x)) with G G Z[A, T] he a nonzero 
Pfaffian function. For every a G M such that g{(y) = 0, we have 

mult(a,ff) < 2 degx(G)degy(G) +degy(G)((5y - 1 ) + {5x + l)degy(G). 
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Proof. Assume first that G is irreducible in 7j[X,Y]. If g{a) = 0, then mult(a, 5 ) > mult(a, 5 '). 

As g'{x) = then G does not divide G and, therefore, R := Resy(G, G) / 0. Let 

S,T £ X[X,Y] be such that R = SG + TG. We have that 

R{x) = S{x,ip{x)).g{x) +T{x,(p{x)).g'{x). 

If a is a multiple root of g, the previous identity implies that mult(a, 5 ) < mult(a, i?) + 

1 < d^(i?) + 1. Taking into accourh that deg(i?) < degjj^(G) degy (G) + degj^ (G) degy (G), 
degjsf(G) < degjsf(G) + 6x and degy(G) < degy(G) — 1 + (5y, we conclude that 

mult(a,fl') < 2degx(G) degy(G) +degx(G)((5y - 1) + 5jvdegy(G) + 1. 

In the general case, write G = c{X) ni<i<t Gi{X, where c{X) = cont(G) and Gi,... ,Gt € 
Z[X,Y] are irreducible polynomials. For every i, let gi{x) = Gi{x,ip{x)). From the previous 
bound, we deduce 

mult(Q:, 5 ') = mult(a, c) + rrii gi) < 

l<i<t 

— degj\:(c) + rrii {2degx{Gi) degy(Gj) + degj^(Gj)((5y — 1) + 6x degy(Gj) + 1) 

l<i<t 

< 2degj^(G) degy(G) + degjc(G)((5y - 1) + {5x + l)degy(G). 

□ 

The theoretical results in the previous section enable us to construct the following algorithm 
for zero counting for a function f[x) = F{x,ip{x)), where F G Z[X,Y]. By Lemma [HI we will 
assume that Resy(F, F) ^ 0. 


Algorithm ZeroCounting 

INPUT; A function (/? satisfying a differential equation (p'{x) = <b(x, (/^(x)), a polynomial F G 
Z[X,Y] such that Resy(F, F) ^ 0, and a closed interval [a,/3] C Dom(99) with a,/3 G Q. 
OUTPUT: The number of zeros of f{x) = F{x,(p{x)) in [a,/3]. 


1. LetFi(A,y) := 


F{X,Y) 


if degy(F) < degy(F) 


Remainder (lc(F)^F, F) otherwise 

est even integer greater than or equal to 1 + degy (F) — degy(F). 


where D is the small- 


2 . Compute the polynomials Ri and Tj, for 0 < z < A, and pj, for 2 < z < A + 1, associated 
to F and Fi as in Notation O 


3. Determine and order all the real roots ai < a 2 < ■ ■ ■ < cuk lying in the interval (a, b) of 
the polynomials Tj, for 0 < i < N, and pi, for 2 < z < A + 1. 

4. For every 0 < j < k, compute the Sturm sequence fj. = {fij,i)o<i<N for / in Ij = (aj, Oj+i) 
as in Definition dUl where ao = a and a^+i = /3- 

5. Decide whether /(oj) = 0 for every 0 < j < /c + 1 and count the number of zeros. 
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6 . For every 0 < j < k, compute Vj := v{ii.,a'^) — 

k 

7. Compute ^^{0 < j < k + 1 : f{aj) = 0} + ^ Vj. 

i=i 


Complexity analysis: 

Let dx '■= degjf(F), dy '■= degy(F) and, as before, 6x ■= degj^($), dy := degy($). 

Step 1. Note that degy(Fi) < dy. In the case when degy(F) > dy., in order to bound 
degjf(-Fi), notice that degjf (lc(-F)^F) < Ddeg(lc(F)) + dx + dx- Then, the polynomial 
Fi can be obtained by means of at most D successive steps, each consisting of subtracting 
a multiple of F with degree in X bounded by {D — i) deg^(lc(F)) + (i + l)dx + dx from a 
polynomial whose degree in X is bounded by (D — i + 1) degji(^(lc(F)) + idx + dx- Then, 
degxCTi) < {D + l)dx + dx < {dy + 2)dx + dx- 

In order to perform the computations (as polynomials in the variable Y) avoiding division 
of coefficients (which are polynomials in X), we do not expand the product of the coeffi¬ 
cients of F times lc(T)'^ at the beginning, and at the ith step, we write each coefficient of 
the remainder as a multiple of lc(-F)^“*. Thus, at each step, we compute at most dy dy 
polynomials in X: for the first dy of them, we compute the difference of two products 
of a coefficient of F (whose degree is at most dx) by a polynomial of degree bounded by 
{i + l)dx + dx, and for the other ones, the product of the leading coefficient of T by a 
polynomial of degree bounded by {i + l)dx + dx- Then, the overall complexity of this step 
is 0{{dy + dy)dxdy{dydx + dx))- 

Step 2. Each subresultant of F and F^ is a polynomial in the variable V whose coefficients are 
polynomials of degree bounded by (dy — l)dx + dy{{dy + 2)dx + dx) in the variable X. 
We compute it by means of interpolation: for sufficiently many interpolation points, we 
evaluate the coefficients of F and Fi, we compute the corresponding determinant (which 
is a polynomial in Y with constant coefficients) and, finally we interpolate to obtain each 
coefficient. 

For each interpolation point, the evaluation of the coefficients of F and Fi can be performed 
within complexity 0{dydx + (dy — 1)(((5y + 2)dx + <5x)) = 0{dy{dydx +5_x))- Then, we 
compute at most 2dy — 1 determinants of matrices of size bounded by 2dy — 2 within com¬ 
plexity 0 (dy+^), we multiply them by the polynomials Y^ F or Y^ Fi evaluated at the point 
and we add the results in order to obtain the specialization of the subresultant at the point, 
which does not modify the complexity order. This is repeated for dy{{dy + 3)dx + dx) 
points. Finally, each of the at most dy coefficients of the subresultant polynomial is com¬ 
puted by interpolation from the results obtained. Each polynomial interpolation can be 
done within complexity 0{M{dy{dydx + dx)) ^og{dy{dydx + dx)))- Then, the computa¬ 
tion of the at most dy coefficients of each subresultant can be achieved within complexity 
0{{dy{dydx + dx) + dY'^^)dy{dydx + dx) + dyM{dy{dydx + dx)) log((iy((Iy(ix + (5x))) = 
0{d‘^+\dydx+dx)‘^)- 

As we have to compute at most dy subresultants, the overall complexity of the computation 
of all the required subresultants is of order 0{dY^^{dydx + dx)"^)- 
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Note that we may compute successively only the polynomials Ri = SReSnj-i- The index 
nj_|_i indicating the next subresultant to be computed is the degree of Ri, and the poly¬ 
nomial Tj is its leading coefficient. Finally, the polynomials pi € are subresultant 

coefficients of F and Fi, which are also computed by interpolation. The complexity of 
these computations does not modify the order of the overall complexity of this step. 

Step 3. Consider the polynomial 


L{x)= n n p^- 

0<i<N ‘i<i<N+l 

Note that P2 = (-l)l(degy(F)-degy(Fi))(degy(F)-degy(Fi)+l)jp(^^)degy(F)-degy(Fi). 
the same zeros as ti = lc(Ti). 

We determine the Thom encodings of the roots of L in the interval (a, h) by computing the 
realizable sign conditions on Der(L), X — a, {3 — X, where Der(L) = (L, L ',..., 

The degree of L is bounded by {2d\ — dY){{6Y + 3)dx + dx)- We compute its co¬ 
efficients by interpolation: the specialization of L at a point can be computed within 
0{dY{dYdx + operations by specializing its factors and multiplying, and this is done 
for deg(L) -|-1 points; then, the total complexity of evaluation and interpolation is of order 
0(dy((5ydx -|- dx)'^)- The complexity of computing the realizable sign conditions is of 
order 0{d^{dYdx + dxY log^((iy (Jydx + dx))) (see Lemma ED . Finally, we can order the 
roots of L in {a, (3) by comparing their Thom encodings (see [ll Proposition 2.28]) within 
complexity 0((iy (Jydx -|- dxY log((iy ((iydjv -|- dx))) using a sorting algorithm. 

The overall complexity of this step is of order 0{dY{dYdx + log^(dy (5ydx -|- dx)))- 

Step 4. The Sturm sequences i^ij)o<j<k are obtained by multiplying the polynomials {Ri)o<i<N 
by the corresponding signs (cr/^-,i)o<i<Af as stated in Definition [TOl Note that if p is a 
univariate polynomial having a constant sign in Ij = (ccj, aj+i), to determine this sign it 
suffices to determine sg(p, a^) or sg{p,ajYi), which can be obtained from the signs of p 
and its successive derivatives at aj or Oj+i respectively. 

Then, in order to compute the required signs, we compute the realizable sign conditions 
on the family 

Der(L), X -a,P - X, DeT{pi)s<i<N, Der(ri)i<i<Ar_i 

which consists of 0{d\{dYdx +<^x)) polynomials of degrees bounded by (2(iy — (iy)(((5y -|- 
3)dx+dx)- The complexity of this computation is of order 0(dy ((5ydx+<^x)^ log^(dy ((5ydy:+ 
dx)))- Going through the list of realizable sign conditions, we determine the signs aj.^i 
and, from them, the Sturm sequences ij. within the same complexity order. 

The overall complexity of Steps 1 - 4 is of order 0{dY{dYdx + log^(dy (5ydjv + (^jv)))- 

Steps 5 and 6. These steps require the determination of the sign of Pfaffian functions of the 
type G{x, ip{x)), with G € ^[X, Y], at real algebraic numbers given by their Thom encod¬ 
ings (more precisely, at the real roots Uj of L lying on (a, /3) and at the endpoints a and 
(3 of the given interval). We assume an oracle is given to achieve this task. 
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At Step 5, we need k + 2 < deg(L) + 2 = 0{dY{5Ydx + dx)) calls to the oracle for 
the Pfaffian function dehned by the polynomial F, having degrees degjj('(-F) = dx and 
degy(F) = dy. 

At Step 6 , we use the oracle for Pfaffian functions defined by polynomials with degrees in X 
bounded by dyddy FS'jdx -hdx) and degrees in Y bounded by dy. Taking into account the 
bound for the multiplicity of a zero of such a function given by Lemma fTSl it follows that the 
determination of sg(//^.^j, ) and otj) requires at most 0{dy{dy+5y){5ydx+dx)) 

calls to the oracle. Then, the oracle is used at most 0{dY[dy + 6y){6ydx + dxY) times. 

Therefore, we have the following: 

Proposition 16 Let f{x) = F{x,(p{x)) be defined from a polynomial F G Z[A, T] and a Pfaf¬ 
fian function y? satisfying ip'{x) = <h(x, (^(x)), where <f> G Z[X, T] with degy(‘h) > 0. Let 
dx '■= degx{F), dy := degy(T), 5x '■= degj(^(<l>) and 6y := degy(<b). Then, Algorithm 
ZeroCounting computes the number of zeros of f in a closed interval [q;,/3] C Dom{ip) (a, fi G Q) 
within 0 ((iy ((5ydx+ <5x)^ log^((iy + <5x))) arithmetic operations and comparisons, and us¬ 

ing at most 0{dd{dy+ dy){6ydx + dx)‘^) calls to an oracle for determining the signs of Pfaffian 
functions of the type G{x,ip{x)), with G G Z[A, T], at real algebraic numbers. 

As a consequence of the previous algorithm we deduce an upper bound for the number of 
zeros of the Pfaffian functions under consideration in a bounded interval: 

Corollary 17 Let f{x) = F{x,ip{x)) be defined from a polynomial F G Z[X, T] and a Pfaffian 
function ip satisfying p>'{x) = <l>(x, where $ G Z[A, T] with degy(<h) > 0. Let dx '■= 

degx(T); dy := degy(F), 6x ■= degx(^) and 5y := degy(<l>). Then, for any open interval 
I C Dom((/?), the number of zeros of f in I is at most {dy + l)(2dy — dy){{5y + Yjdx + 5x)- 

An alternative bound can be obtained from Khovanskii’s upper bounds for the number of 
non-degenerate zeros of univariate Pfaffian functions and for the multiplicity of an arbitrary 
zero of these functions (see [5]). Keeping our previous notation, for a polynomial F G ’L\X,Y] 
with deg(T) = d, if deg(<k) = 6, using Khovanskii’s bounds, it follows that both the number of 
non-degenerate zeros and the multiplicity of an arbitrary zero of f{x) = F{x,ip{x)) are at most 
d{5 + d). We can get an upper bound for the total number of zeros of / by bounding the number 
of non-degenerate zeros of / and of its successive derivatives of order at most d{6 + d) — 1. 

Following dH), we have that f is dehned by a polynomial of degree at most d + 6 — 1 and so, 
for every /c G N, is given by a polynomial of degree at most d -|- k{5 — 1). Then, the total 
number of zeros of / is at most 
d{5+d)-l 

{d + k{6 -l)){d + d + k{d - 1 )) < ^d^6^{d + df. 
k=o ^ 

Note that the bound from Corollary [T3 is of lower order than this one. 

5 E-polynomials 

In this section, we will deal with the particular case of K-polynomials, namely when (p{x) = 
for a polynomial h G Z[A] of positive degree. We will hrst show how to perform steps 5 and 6 of 
Algorithm ZeroCounting (that is, we will give an algorithmic procedure to replace the calls to 
an oracle). Finally, we will prove a bound for the absolute value of the zeros of an Fl-polynomial. 
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5.1 Sign determination 

The main goal of this section is to design a symbolic algorithm which determines the sign that 
an ill-polynomial takes at a real algebraic number given by its Thom encoding. To do this, we 
will use two subroutines. The first one, which follows m Lemma 15], determines the sign of an 
expression of the form — a for real algebraic numbers a and /3. The second one allows us to 
locate a real number of the form for a real algebraic number a, between two consecutive 

real roots of a given polynomial. 


Algorithm SignExpAlg 

INPUT; Real algebraic numbers a and /? given by their Thom encodings ap^{a) and ap2{j3) 
with respect to polynomials Pi,P 2 € Z[X] such that deg(Pi), deg(P 2 ) < d (d > 2) and 
H{Pi),H{P2)<H. 

OUTPUT: The sign s := sign(el^ — a). 

1. Let c := l)iL)-2"««®(5««+4riog(H)l)_ 

2 . Compute w ^ Q such that je^ — w\ < c as follows: 

(a) Compute rci G Q such that \/3 — wi 

(b) Compute w £ Q such that — w 

3. Compute s = sign(r(; — a). 


2. 3^+2 


c 



Proof of correctness and complexity analysis: 

Step 1. We will show that, for the chosen value of c, the inequality |e^ — a| > c holds. 

As shown in m, if a and (3 are algebraic numbers of degrees bounded by 6 and heights 
bounded by v, then 

- aj > e-2'‘^'^®l°(^+e0(ln(0+lnln(i/)) 

Note that 

g 2 « 06 ln(i.+e-)(lnM+lnlnH) < ^ < (j, + 16)2^®®® 

It is clear that the degree of an algebraic number is bounded by the degree of any polyno¬ 
mial which vanishes at that number. With respect to the height, by [U Propositions 10.8 
and 10.9], we have 

H{a) < 2‘^||PiII < 2^{d + 

and, similarly, it follows that the same bound holds for Here, ||Hi|| stands for the 

norm 2 of the vector of the coefficients of Pi. 

The required inequality is deduced by taking 6 = d, u = 2‘^(d -|- and using the 

bounds 

2 '^(d + l)^/ 2 ^ + 16<2'^+^(d+l)H and ln( 2 '^(d + l)^/^^) < ^d + [log(P)l. 
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Step 2.(a) Applying the algorithm from Lemma S] to the polynomial P 2 with e = „ , we 

3a+d 

get intervals Ij = {aj,hj\ with aj,bj G Q and bj — aj < e {1 < j < k) such that /3 G Ij^ 
for some jo- We determine the index jo by computing the feasible sign conditions for 
Der(P 2 )j A1 — ai, X — bi,..., X — X — 6 ^- Finally, we take wi = bj^. The complexity of 
this step is of order .c~^) + log^(d))) = 0{(PH + dP{d + log(Lf))^). 


By the mean value theorem, the inequality |/3 — rcil < 


2. 3^+2 


implies that |e^ — 


c 

< 2 - 


Step 2.(b) Following iini Lemma 14], in order to obtain w, we compute the Taylor polynomial 
centered at 0 of the function e® of order t := 8([log(2/c)] + 1 + H) specialized in wi. The 
complexity of this step is bounded by 0{(f{d + log{H))‘^ + H). 


Step 3. The fact that sign(u) —a) = sign(e^ —a) is a consequence of the inequalities \e^ — a\ > c 
and |e^ — 'u;| < c. In order to determine this sign, we compute the feasible sign conditions 
on Der(Pi), A — w and look for the one which corresponds to the Thom encoding of a. 
The complexity of this step is of order 0{d^ log^(d)). 


The overall complexity of this subroutine is 0{d^H + dP{d + log(Lf))2). 


The second subroutine is the following: 


Algorithm RootBox 

INPUT: A polynomial h G Z[A], an algebraic number a G M such that h{a) 7 ^ 0, given by its 
Thom encoding as a root of a polynomial L G Z[A], and a polynomial M G Z[A] together with 
the ordered list of Thom encodings of all its real roots Ai < A 2 < • • • < Am- 

OUTPUT: The index io) 0 < *0 < such that A^g < 6^^“^ < Ajg+i, where Aq = — 00 and 
Am+l = + 00 . 

1. Compute S{T) := Resx{L{X),T - h(X)). 

2. Compute the feasible sign conditions on Der(L), S{h), S'(h),..., and the Thom 

encoding of h{a) as a root of S. 

3. Compute sign(e^(“^ — A*) applying Algorithm SignExpAlg, for i = 1,..., m, until the hrst 
negative sign is obtained for zq. If all the signs are positive, zq = m. 


Proof of correctness and complexity analysis: 

Note that h{a) is a root of the polynomial S G Z[T] computed in Step 1. Therefore, in Step 
2, the sign condition on Der(L), S{h), S'{h ),..., (h) having the Thom encoding of a as 

a root of L in the hrst coordinates has the Thom encoding of h{a) as a root of S in the last 
ones. 

Assume that deg(L) < £, deg{h) < 6 and deg(M) < rj. 
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The resultant computation in Step 1 can be done within complexity by interpola¬ 

tion, noticing that deg(5) < 1. Applying LemmaO the complexity of Step 2 is 0{i^6 log(^) log^{£d)). 
Finally, taking into account that H{S) < {£ + 6)1 H{LY{2H{h)Y, defining 

n := max{H{M), {£ + 5)! H{lY i2H{h))Y, 

the complexity of Step 3 is o(jnmaK{r], £Y (ji -|-maxjry, £}®(max{? 7 , £} -|-log(?^))^^y 

The overall complexity of the algorithm is of the same order as the complexity of Step 3. 

Now we are ready to introduce the main algorithm of this section. 

Algorithm E-SignDetermination 

INPUT: Polynomials G G Z[X,V], h G Z[A'], deg(h) > 0, L G Z[A] and Thom encodings 
) <^L(o!t) of real roots oi,..., of L. 

OUTPUT: The signs of G(aj, for l<j<t. 

1. For every 1 < j < t, determine whether G{aj,Y) = 0. If this is the case, the sign of 

is 0. 

2. Compute R = gcd(L,/i) and the list of realizable sign conditions on Der(L), i?, G(X, I). 
Going through the list, determine the sign of G{aj,e^^°‘^Y — G{aj, 1) for every j such that 
G{aj,Y) ^ 0 and R{aj) = 0. 

3. Compute M{Y) := Resx{L{X),G{X,Y)). 

4. Compute the Thom encodings of the real roots of M and order them: Ai < • • • < Am- 

5. For every 1 < j < t such that G{aj, Y) ^ 0 and R{aj) / 0: 

(a) Determine the index 0 < ij < m such that Xi^ < < Ai^.+i by applying subroutine 

RootBox, where Aq := —oo and Am+i := -|-oo. 

(b) Find wj gQD (A,^., Ai^.+i). 

(c) Compute the sign of the polynomial G{X,Wj) at X = aj. This is the sign of 


Proof of correctness and complexity analysis: 

Assume that degji(^(G) < dx, degy(G) < dy, deg(L) < £ and deg(h) < 6. 

Due to Lindemann’s theorem, if a € M is an algebraic number and h{a) 0, then is 
transcendental over Q. Therefore, for an algebraic number a G M, G(a, = 0 if and only if 

either G{a,Y) = 0 or h{a) = 0 and G{a, 1) = 0. Then, Steps I and 2 enable us to determine 
all the indices j such that G{aj, = 0. 

Step 1. Compute cont(G), the gcd of the coefficients of G, by applying successively the fast 
Euclidean algorithm [3l Algorithm 11.4] within complexity 0(dyM(dx) log(dx))- Then, 
determine the realizable sign conditions on Der(L), cont(G) within max{£, dx} log(^) 
log^(max{I', dx})) arithmetic operations. 
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Step 2. The complexity of the computation of R is of order 0(M(max{^, 5}) log(max{£, <5})) 
and the realizable sign conditions on Der(L), i?, G(X, 1) can be found within complexity 
0{£‘^ max{^, dx}^og{t) log^(max{£, dx}))- 

Step 3. In order to compute M{Y), evaluate G{X,y) at sufficiently many values y, compute 
the corresponding determinants and interpolate. Taking into account that deg(M) < idy, 
the total cost of this step is of order 0{idY{dx + + M{idY) log{£dY))- 

Step 4. The computation of the Thom encodings of the real roots of M can be done within 
0{{^dYY^o^{^dY)) operations. Then, we order the real roots of M by means of their 
Thom encodings within complexity of order 0((£dy)^ log(.£(iy)). 

Step 5. Following the proof of [U Proposition 8.15], it follows that H{M) < {I + dx^iidy + 

Recall that deg(M) < Idy- 

(a) The complexity of this step is 0{{£dyY{'H + {idy^i^dy + log{'H)Y)), where Ti = 
max{(^ + 6)\H{Lf{2H{h))Y {£ + dx)\H{LY^ {{dy + l)H{G)f}- 

(b) By applying Lemma 0] to the polynomial M and a lower bound e for the minimum 
distance between two different roots of M, we obtain pairwise disjoint intervals (a*, 6 j] 
with rational endpoints such that A* G (oi, 6 i] for i = 1,..., m. Following Lemma El 

we can take e = (idy)-{idy + 1) {{i + dx)! H{LY^ {{dy + l)H{G)Yy-^^^ . 

Let Wj := bi-. 

The complexity of this step is 0((-^dy)^((^+dx) log(-^+dx)+^(log(LI(G))+log(dy)) + 
dxlog(id(L)))). 

(c) We compute the coefficients of G(X,Wj} within complexity 0{dxdy). Then, we com¬ 
pute the feasible sign conditions of Deic{L),G{X,Wj), which enable us to determine 
the sign of G{aj,Wj), within 0{i'^ max{^, dx}log{i) log^(max{^, dx}))) additional op¬ 
erations. 

The overall complexity of the algorithm is 0{t{£dYY{'H + {idyYYdy -|- log(^))^)). 

The previous complexity analysis leads to: 

Proposition 18 Given polynomials G G Z[X, T], h G Z[X], deg(/i) > 0, L G Z[X] with 
degrees bounded by d and height bounded by H, and Thom encodings ^^(ai),..., aL^ott) of real 
roots of L, we can determine ^{1 < j < t : G{aj,e^^°^^Y — 0} within complexity 

0(d^ log^(d)). Moreover, the signs of G{aj, for 1 < j < t, can be computed within 

complexity 0{t8'^d^'^~^^H‘^Y- 

5.2 Zero counting for E'-polynomials 

Here, we will apply Algorithm E-SignDetermination from the previous section as a subroutine 
in Algorithm ZeroCounting described in Section 0] to obtain a zero counting algorithm for 
Fi-polynomials with no calls to oracles. 

In order to estimate complexities we will need upper bounds for the degrees and heights of 
polynomials defining the successive derivatives of an Fi-polynomial. 
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Remark 19 For a Pfaffian function g{x) = given by a polynomial G G Z[X, y], we 

~ ~ dG dG 

have that g'{x) = G{x,e^^^^) is given by the polynomial G := + h'{X)Y-^. If degxiG) = 

dx, degy(G) = dy and deg(/i) = 5, we have that 

degx(^) < 5 — 1 + dx, degy(G) = dy 
H{G) < H{G){dx + dy6^H{h)) 

Applying these bounds recursively, we get that the successive derivatives of g can be obtained as 

for polynomials G Z[X,Y] sueh that 

degxC'G) < u{5 -l)+dx, degy{''G) = dy 
u-l 

HC^G) < H{G) -l)+dx + dy6^H{h)). 

j=0 

Now, we can state the main result of this section. 

Theorem 20 Let f{x) = F{x,e^^^'l) be an E-polynomial defined by F ^ X[X,Y] and h G Z[y] 
with deg(F), deg(/i) < d and H{F), H{h) < H, and let [a, 6] he a closed interval. Assume that 
Kesy{F,F) 7 ^ 0. There is an algorithm that computes the number of zeros of f in [a,b] within 
complexity 

Proof. In order to prove the theorem, we adapt Algorithm ZeroCounting introduced in Section 
mto count the number of zeros of an ill-polynomial with no call to oracles. It suffices to show 
how to perform Steps 5 and 6 of the algorithm and estimate the complexity of the procedure. 

Step 5 can be achieved by means of Steps 1 and 2 of Algorithm E-SignDetermination. As 
in this case deg(L) < lOd^, the complexity is of order 0{d^ log^{d)). 

To achieve Step 6 of the algorithm, we apply the algorithm E-SignDetermination to the 
polynomials defining the functions fj.^i and their successive derivatives, for 0 < i < N. These 
functions are defined, up to signs, by the polynomials Ri introduced in Notation [9l and ^Ri, 
0 < i < N, G N. 

Since degy(F) = degy(F), then Fi = lc{F)‘^.F — lc{F)lc{F)F and so, degy'(Fi) < Ad—l and 
H{Fi) < 4:d{d-\- l)H^{d-\-d^H) < 8((i-|- l)d'^H'^. Taking into account the determinantal formula 
for the subresultants, it follows that for every k, degjj^ (SResfc) < 5d^ — 2d and R(SResfc) < {2d — 
l)!25(i-2(^- 1 - -^'^2d-2^5d-i j^5d-i ^ ^2d-i25d-2^9d-3 jjbd-i^ which are therefore, upper bounds for 

degy(i?j) and H{Ri) for all i. Finally, recalling that L is the product of at most 2d polynomials 

of degrees at most 5d^ — 2d that are coefficients of the subresultants SRes^, we have that 
H{L) < (^I,(f^‘^d-1^^2d-l25d-2^9d-3jj5d-l\^2d ^ ^4d^ 2 ^^'^'^^18d'^-2d-2^lOd^-2d 

Taking into account the bound for the multiplicity of a zero of a Pfaffian function from Lemma 
da we will apply the algorithm E-SignDetermination to the polynomials Ri {0 < i < N) and 
'^Ri for V < lOd^ — 3d^, to determine the signs of the corresponding Pfaffian functions at the zeros 
of L. The bounds from Remark 1 19 1 applied to the polynomials Ri imply that, for v < lOd^ — 

degxi''Ri) < (lOd^ - 2,d^){d - 1) + hd^ -2d< lOd^ - 5d^ 
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HC'Ri) < H{Ri){Wd'^ + iH - 

Then, the complexity of applying the algorithm to each of these polynomials is of order 

Oid^^in + d^^id"^+ log'Hf)) 

where 

H < (lOd^ + 5d^)lH{L)^^^*-^^\{d + + {H - 

= {2dH)^^'^^\ 

This sign computation is done for at most d{lQd^ — 3d^) polynomials. Finally, for each 
interval Ij, the signs sg{fi.^i,a^) and sg{fi.^i,ajj^i) are obtained easily following Definition fTOl 
Therefore, the overall complexity of the algorithm is of order 

{2dH)°^'^^\ 


□ 


The previous procedure can be slightly modified to count algorithmically the total number 
of real zeros of an i?-polynomial. To do this, we consider the signs of Fi-polynomials at +oo and 

—CX). 

Let g{x) = G{x, be an £'-polynomial. Assume G{X, Y) = J2j=o with Udy / 0 

and let jo = min{j : aj ^ 0}. We define 


sg(5, Too) 


sign(lc(ajp)) if lc(/i) < 0 
sign(lc(ady)) if lc(/i) > 0 


and 

_ fsign((-l)‘^®g(“™)lc(ajJ) if (-l)d'^gWlc(/i) < 0 
Sg(ff, -oo) - > q 

For a sequence of Fi-polynomials f = {fo,..., Jn), we write u(f,+oo) for the number of 
variations in sign in (sg(/o, +oo),..., sg(/Ar, +oo)) and u(f, —oo) for the number of variations in 
sign in (sg(/o, -oo),..., sg{fN, -oo)). 


Remark 21 Following Notation\E and Deftnition[T(K let and f/_^ be Sturm sequences for 
f{x) = F{x, in the intervals I+oo = {M, +oo) and /_oo = (—oo, —M) where M is an upper 

bound for the absolute values of the roots of Ti for i = 0,..., N and pi for z = 2,..., + 1. 

Then, the number of zeros of f in I+oo equals u(f, M+) — u(f. Too) and the number of zeros 
of f in /_oo equals u(f, —oo) — u(f, —M~). 


By applying this remark, we conclude that the total number of zeros of an S-polynomial in 
M can be determined within the same complexity order as in Theorem 1201 

Remark 22 The assumption ResyiF, F) ^ 0 in Theorem can he removed by using the 
construction in the proof of Lemma 0 Taking into account the increase of height and degree, it 
follows that the overall complexity of the root counting algorithm is of order {2dHY°^^^ as stated 
in TheoremUl 
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5.3 Bound for the size of roots 


The following proposition provides an interval which contains all the zeros of an ii^-polynomial 
and whose endpoints are determined by the degrees and heights of the polynomials involved 
in its definition. Using this bound, applying successively our algorithm for zero counting, it is 
possible to separate and approximate the roots of an £^-polynomial. 


Proposition 23 Let f{x) = be an E-polynomial defined by F ^ Z[X, T] and h G 

Z[X] such that deg(F) < d, deg{h) = 6 > 0 and H{F),H{h) < H. Then, for every zero a € M 
of f, we have that |a| < M{d, 6, H) := 1 + (d + l)id^ max{(d + 1)(1 + 2[^ + Ij!}. 

Proof. Let F{X,Y) = ^ with deg(aj) < dx for every 0 < j < dy and 

Let a G M be a zero of /. If ady^ot) = 0, then |a| < r^ady) < I + H (see Lemma [5]) and so, 
the bound in the statement holds. Similarly, if ao(a) = 0, the bound holds. 

Assume now that arfy(a) 0 and ao(a) 0. Then is a root of F{a,Y) G M[y] and 
g-h(a) -g q£ F{a,Y~^) G M[y]. By Lemma [21 it follows that 



We are going to prove that, for a > M{d, S, H), one of the previous inequalities fails to hold. 
Note that in both cases, the right hand side of the inequality is given by a rational function. 


E 


0<j<dy 


aj{Xf 


adyiXfi 


and 


E 


0<j<dy 


aj{Xy^ 


aoixy 


respectively, where the numerator and the denominator are integer polynomials of degrees at 
most 2dx and coefficients of size bounded by (dy + l)(dx + l)dd(F)^ and (dx + l)dd(T)^ 
respectively. Moreover, the degree of the denominator is less than or equal to the degree of the 
numerator. 

First, assume that the leading coefficient of h is positive. 

Let p{X) = Eo<j<dy and q{X) = aj^iX) so that ^ = 1 + Eo<i<dy-i • 

and let C > 0 be the quotient of the leading coefficients of p and q. Note that [Cl < (dy + 


l)H{Fy. 

If deg(p) = deg(g), for every x > max{r{q),r{p — (C + !)<?)}, we have that 


Pjx) 

q{x) 


<C + l. 


On the other hand, for x > r{2h — ln(C + 1)), we have that > C + 1. We conclude that, 

for x > m.ax{r{q),r{p — iC + l)q),r{2h — ln(C' + 1)}, the inequality ^ holds. 

q{x) 

If deg(p) > deg(g), let do := deg(p) — deg(g). For x > max{r(g), r(p —2C'x'^°g)}, we have that 
< 2Cx^°. Note that for x > r{2h — X^). As e^^ > ^ 2Cx^° 

Q[X ) rC! 

for x > ^’(Ei.-n'''^^ — 2CX'^°), it follows that < g 2 h{x) ^ ^ max{r(g),r(p — 

^ k\ q{x) 

2Cx'^°q),r{Ylkl(^^^ ~ 2C'A'^°)}. Using again Lemma[2l we obtain; 
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. r{q) <l + {dx + 1)F(F)2 

. rip -iC + l)q)) <l + idx + l)HiF)^dY + [dy + l)HiFf) 

• r(2/i - ln(C + 1)) < 1 + Hih) + ^ ln((dr + l)ff(F)^ + 1) 

• r(p - 2CX^°q) <l + idx + l)(dy + l)HiF)^l + 2HiF)^) 

• ri2h-X^) < l + 2Hih) 

• - 2C7X'^o) <l + 2Y^ + l\\idY+ l)H{Ff 

and, therefore, we conclude that, for a > M(d, <5, if), the following inequality holds 



If the leading coefficient of h is negative, applying the previous argument to —/i, we have 
that, for a > M{d,6,H), the following inequality holds 



Finally, noticing that a is a zero of if and only if —a is a zero of F(—x,e^*' 

we conclude that every zero a of f satisfies a > —Mid,5,H). □ 
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